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The X-ray afterglow flat segment in short GRB 05 1221 A: 
Energy injection from a millisecond magnetar? 
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ABSTRACT 

The flat segment lasting ~ 10 4 seconds in the X-ray afterglow of GRB051221A rep- 
resents the first clear case of strong energy injection in the external shock of a short 
GRB afterglow. In this work, we show that a millisecond pulsar with dipole magnetic 
field ~ 10 14 Gauss could well account for that energy injection. The good quality 
X-ray flat segment thus suggests that the central engine of this short burst may be a 
millisecond magnetar. 

Key words: Gamma Rays: bursts-GRBs: individual (GRB 051221 A) -ISM: jets 
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1 INTRODUCTION 



GRB 051221A was localized by the Burst Alert Telescope 
(BAT) onboard the Swift satellite (Parsons et al. 2005) and 
promptly observed by both Swift /BAT and the Konus-Wind 
instrument. The Swift observations reveal this is a short 
hard burst, with T90 = 1.4 ± 0.2 s, a hard photon index 
q = -1.39±0.06, and a fleunce 1.16±0.04 x 10" 6 ergscm~ 2 
in the 15-150 KeV band (Cummings et al. 2005). The Konus- 
Wind cutoff power-law spectral fitting, in the 20-2000 KeV 
band, shows a fluence 3.2^1' 7 X 10 -6 erg cm -2 , a low-energy 
photon index a = —1.08 ± 0.14, and an observed peak en- 
ergy S pcak = 402±?2 KeV (Golenetskii et al. 2005). With 
a redshift z = 0.5459 (Soderberg et al. 2006), this burst's 
isotropic prompt emission energy is _E 7 ~ 2.4 x 10 51 erg, us- 



0.27, n 



A 



0.73, 



ing the ACDM concordance model of Am 
and h = 0.71. 

Both the X-ray (~ 10 2 - 2 x 10 6 s) and the optical 
(~ 10 4 - 4 x 10 5 s) afterglow light curves of GRB051221A 
have been well detected, while in the radio band only one 
data point followed by several upper limits is available. This 
burst is distinguished by an X-ray flattening at t ~ 0.03 — 0.2 
day, which strongly suggests a significant energy injection 
(Soderberg et al. 2006; Burrows et al. 2006). However, the 
nature of that energy injection is not clear. In the widely ac- 
cepted double neutron star merger model for the short/hard 
burst, supported by the lack of detection of the bright super- 
nova component in the current event (Soderberg et al. 2006), 



the material ejected in the merger is ~ (10 -4 — 1O~ 2 )M0 
(Rosswog et al. 1999; Ruffert & Janka 2001). Given an en- 
ergy conversion efficiency ~ 0.001, the fall-back accretion 
of part of that material onto the central post-merger ob- 
ject is not likely to be able to pump energy up to ~ 10 52 
erg, even with moderate beaming correction. So the fall- 
back accretion model, which may give rise to significant en- 
ergy injection in the collapsar scenario of long/soft GRBs 
(MacFadyen, Woosley & Heger 2001), does not work in the 
current case. 

In this Letter, we'll show that the afterglow undergoing 
an energy injection from a millisecond pulsar with a dipole 
magnetic field ~ 10 14 Gauss (i.e., a magnetar) well accounts 
for the multi-wavelength data. The good quality X-ray flat 
segment thus suggests that the central engine of this short 
burst may be a millisecond magnetar. 



2 ANALYTICAL INVESTIGATION 

In the long/soft GRB scenario, the energy injection caused 
by a millisecond pulsar has been discussed in some detail 
(Dai & Lu 1998; Wang & Dai 2001; Zhang & Meszaros 2001; 
Dai 2004; Ramirez-Ruiz 2004; Zhang et al. 2006). Similarly, 
provided that the gravitational wave radiation is not impor- 
tant, the dipole radiation luminosity of a magnetar can be 
estimated by 



L dip (t b ) ~ 2.6 x 10 48 erg s" 1 B\ 



4 < 6 ^(l+^/T )- 2 ,(l) 
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where B± is the dipole magnetic field strength of the magne- 
tar, R s is the radius of the magnetar, f2 is the initial angular 
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frequency of radiation, the subscript "6" represents the time 
measured in the burst frame, T = 1.6 x 10 4 B^ 2 14 O4 2 /45_R~g 
s is the initial spin-down timescale of the magnetar, / ~ 
10 45 g cm 2 is the typical moment of inertia of the magnetar 
(Pacini 1967; Gunn & Ostrikcr 1969). Here and throughout 
this text, the convention Q x — Q/10 x has been adopted in 
cgs units. 

The energy emitted at t b will be injected into the pre- 
vious GRB ejecta at time T b satisfying 



[1 - (3(T b )]cdT b = C / (3{T b )dT b » Ctb 



tb 



I 

Jo 



(2) 



where (i, in units of the speed of light c, is the velocity of the 
ejecta moving toward us. The corresponding observer time 
is 



stitute p and / into equation (7), we find Ek ~ 7.6 x 
10 745 ^4 er g is we ll consistent with the observational result 
E k ~ 2E-, ~ 4.8 x 10 51 erg when taking typical I45 ~ 1.5, 
Q4 ~ 0.65 (i.e. 1 millisecond period), and the GRB effi- 
ciency r\ — Ek/(Ek + E 7 ) ~ 30%. Furthermore, the mea- 
surements of t c , tf, /, and p are well consistent with the 
constraint relation / ~ (t//t c )^ p+2 ^ 4 - So we conclude that 
the central engine may be a magnetar. Its physical parame- 
ters are (Q,R 3 ,B ± ,I) ~ (6500 s"\ 13 km, 10 14 Gauss, 1.5 x 
10 45 g cm 2 ) according to the above constraint relations. 

An additional constraint is on the ellipticity e of the 
magnetar. As shown in Shapiro & Teukolsky (1983), the 
spin-down timescale due to the gravitational wave radiation 
is r gw - 3 x 10" 3 e" 2 /^ 5 1 n4 4 s. Now r gw > t f ~ 1.5 x 10 4 s, 
we have e < 5 x lO" 4 /^ 172 ^ 2 . 



t« (1 + . 



[1 - p(n)]dn. 



(3) 



Equations (2) and (3) yield + z)t b + (1 + z) J*"[l - 

f3(Tb)]dTb « (l + z)t b . 

At time T b , the energy injected into the ejecta satisfies 
dE/dT b = [1 - /3(T b )]L dip (t b ). With the relation dt = (1 + 
z)[\ - (3(T b )]dT b , we have 

dE 
dt 



l + z 
2.6 x 10 4? 
(1 + *) 



T + 2 



erg s- 1 B 2 ±M Rl 6 Qi[l + 



(l + z)T ' 



2 (4) 



So the energy injection rate dE/dt ~ const for i <C (1 + «)T 
and dS/dt oc t~ 2 for £ > (1 + z)T a . 

A general energy injection form can be written as 
dE/dt = A(l + z) _1 (t/t jr 9 for U < t < t f , where U and 
t/ are the times when the energy injection takes place and 
turns off, respectively (Zhang & Meszaros 2001; Zhang et 
al. 2006). GRB ejecta's dynamical evolution, at a time t c , 
is significantly changed when the injected energy roughly 
equals to its initial kinetic energy, i.e., J " (dE/dt)dt ~ Ek- 

We thus derive ^(<c" q -*i" q ) ~ {l + z)\l-q)E k . Accord- 
ingly, our magnetar model requires t = 1, U ~ and q ~ 
for t c < (1 + 2)T , which leads to ^4 ~ (1 + z)E 1 /t c , and so 



2.6 x 10 48 erg s" 1 Bl M R%Q.l ~ (1 + z)S fc /t c 



(5) 



For t > (1 + z)T , the rate of the energy injection drops 
sharply or even the central supermassive magnetar has col- 
lapsed when it has lost significant part of the angular mo- 
mentum, which indicates that the afterglow lightcurve flat- 
tening weakens at the time tf > (1 + z)T . So we have 



45 R s , 6 



tf- 



(6) 



1.6 x io 4 (i + 2)b;_ 14 ^^4, 

/.From equations (5) and (6), the total injected energy could 
be estimated as i?mj = tfEk/t c ~ 5 x 10 52 erg Lts^ 2 . On the 
other hand, with and without energy injection, the contrast 
of forward shock X-ray emission flux can be estimated by 
(Kumar & Piran 2000) 



/ ~ (Ei n j/Ek) 



(P+2)/4 



[5 x 10 52 erg U&ljEk 



(p+2)/4 



(7) 



where p ~ 2.4 is the power-law index of the shocked elec- 
trons. Equations (5) — (7) are our main relations to constrain 
the physical parameters of the underlying magnetar. 

From the X-ray observations of GRB051221A, we mea- 
sure t c ~ 3000s, tf - 1.5 x 10 4 s, and / ~ 6. Sub- 



3 NUMERICAL FIT TO THE AFTERGLOWS 
OF GRB 051221A 

We interpret the apparent flattening in the X-ray light curve 
being caused by an energy injection from the central magne- 
tar. Yet the flattening episode is unapparent in the optical 
and radio bands because of limited observations. The code 
here to fit the multi-band lightcurves has been used in Fan 
& Piran (2006) and Zhang et al. (2006), and has been con- 
firmed by J. Dyks independently (Dyks, Zhang & Fan 2006). 
The dynamical evolution of the outflow is calculated using 
the formulae in Huang et al. (2000), which are able to de- 
scribe the dynamical evolution of the outflow for both the 
relativistic and the non-relativistic phases. One modification 
is that now we have taken into account the energy injection 
(for instance see equation (12) of Fan & Piran (2006); see 
also Wei, Yan & Fan 2006). The electron energy distribu- 
tion is calculated by solving the continuity equation with 
the power-law source function Q = K^ p , normalized by a 
local injection rate (Moderski, Sikora & Bulik 2000). The 
cooling of the electrons due to both synchrotron and inverse 
Compton has been considered. The synchrotron radiation 
of the forward shock electrons on the "equal arriving sur- 
face" (on which the emission reaches us at the same time) 
has been calculated strictly. The synchrotron self-absorption 
has also been taken into account strictly. 

We consider a uniform relativistic jet undergoing the 
energy injection from the central source and sweeping up 
its surrounding uniform medium. The energy injection has 
been taken as dE/dt = 2x 10 48 erg s _1 (l+t/1.5 x 10 4 )~ 2 for 
t < 1.5 x 10 4 s otherwise dE/dt = (i.e., we assume that the 
supermassive magnetar collapses when it has lost significant 
part of its angular momentum). As usual, the fractions of 
shock energy given to the electrons and the magnetic field 
(i.e., e e and sb) are assumed to be constant. Shown in Figure 
1 is our numerical fit in the X-ray and optical bands with 
the following jet parameters: Ek = 10 52 erg, e e — 0.3, es = 
0.0002, the circumburst density n = 0.01 cm , the half- 
opening angle 0j = 0.1, and the viewing angle 8 o b s —0 (i.e. 
on-beam viewing). 

In addition, we find that a reverse shock emission 
component, besides the forward shock emission, should be 
evoked to account for the radio data, which is consistent 
with the finding in Soderberg et al. (2006). In the pul- 
sar/magnetar energy injection scenario, the reverse shock 
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Figure 1. Modeling the XRT and R-band afterglow light curves 
of GRB 051221A with energy injection from an underlying mil- 
lisecond magnetar. The optical and X-ray data are taken from 
Soderbcrg ct al. (2006) and Burrows et al. (2006), respectively. 
The r' — band extinction of our Galaxy ~ 0.19 mag and that of 
the host galaxy ~ 0.2 mag (Vcrkhodanov et al. 2000) have been 
taken into account. The solid lines are our numerical results. Fit- 
ted parameters of the GRB jet undergoing energy injection are 
presented in Section 3. 



emission has been calculated in Dai (2004), assuming that 
the outflow is electron/positron pairs dominated and the re- 
verse shock parameters are similar to, or even larger than 
that of the forward shock. The resulting reverse shock emis- 
sion is so bright that long-time flat bump should be evi- 
dent in multi-band afterglow lightcurves. This prediction is 
not confirmed by current observations. This puzzle, however, 
could be resolved if the fraction of the reverse shock energy 
given to the electrons is very small, as shown below. 

In our numerical code, the reverse shock dynam- 
ics/emission has not been taken into account. It is inves- 
tigated analytically instead. We assume that the reverse 
shock emission accounting for the radio data is powered 
by the energy injection. Following Sari & Piran (1999), af- 
ter the reverse shock crosses the ejecta at t x (noted that 
at t ~ t x ~ 0.2 day the central magnetar collapses), the 
observed reverse shock emission flux F ve cx {t/ty.y 2 for 
i/ m < ^® < v c , where v m is the typical synchrotron emis- 
sion frequency of the shocked electrons and v c is the cool- 
ing frequency. On the other hand, _F„ e oc (t/t x ) for 
v® < v-m < v c . The current radio observation thus requires 
a v m (t x ) < 8.46 GHz. At t x ~ 0.2 day, the ejecta is at a 
radius R ~ 7 x 10 17 cm and moves with a bulk Lorentz fac- 
tor r(£ x ) ~ 20 (obtained in our numerical calculation), thus 
the comoving toroidal magnetic field of the magnetar wind is 
B w ~ [2L dip /(-R 2 ric)] 1/2 ~ 2Q/Y W Gauss, where Y w is the 
bulk Lorentz factor of the magnetar wind. In the reverse 
shock phase, the toroidal magnetic field will be amplified 
by a factor ~ T rsh « Y w /2Y{t x ) for Y w > T(t x ) (Kennel & 
Coronitti 1984). On the other hand, the constraint v m (t x ) ~ 
2.8 x 10 6 Hz 7m 2 t ) T(t x )r rah B w /(l + z) < 8.46GHz yields 



7m(i x ) < 20, where 7m(t x ) is the minim random Lorentz fac- 
tor of the electrons accelerated in the reverse shock front. 
Such a small 1 7m(t x ) = 3e e (p - 2)(r rsh - l)m p /[4(p - l)m e ] 
requires that the fraction of shock energy given to the elec- 
trons e e is in order of m e /m p because it is very likely 
that now Y w 3> Y(t x ), where m e and m p are the rest 
mass of the electron and the proton, repsectively. Taking 
e e ~ m e /m p and p ~ 2.3, we thus have r rs h < 115 and 
F w ~ 2r rs hr(t x ) < 4600. After the reverse shock phase, 
the magnetic energy may be translated to the forward shock 
mainly by the magnetic pressure working on the initial GRB 
ejecta and the shocked medium. However, the details are far 
from clear. 

At t ~ 0.91 day, the 8.46 GHz emission flux 
is 0.155 mjy. Therefore the optical (4.6 x 10 14 Hz) 
flux at t x ~ 0.2 day is < 0T55mJy(0. 91/0. 2) 2 (4.6 x 
10 14 Hz/8.46GHz)-°' 65 = 2.7 ^Jy- It is about one 
order lower than the forward shock optical emission. 
As for the X-ray, now v c (t x ) ~ 2 x 10 15 Hz, 
the X-ray (at 2 keV) flux at f x ~ 0.2 day is 
< 0.155mJy(0.91/0.2) 2 (2 x 10 15 Hz/8.46GHz)~°' 65 (4.84 x 
10 17 /2 x lO 15 ) -115 ~ 10~ 6 mjy, is also significantly lower 
than the forward shock X-ray emission. Our estimates made 
here are independent of the poorly known magnetized re- 
verse shock dynamics. We thus conclude that the reverse 
shock X-ray/optical emission is unimportant and can be ig- 
nored. 

One caveat is that the dipole radiation of the magnetar 
is almost isotropic. The medium surrounding the magne- 
tar but out of the GRB ejecta cone would be accelerated 
by the energetic wind. Because the energy injected into the 
GRB ejecta is larger than that contained in the initial ejecta, 
the outflow contributing to the afterglow emission would be 
close to be an isotropic fireball rather than a highly jetted 
ejecta. So the magnetar energy injection model is somewhat 
challenged by the late jet break detected in GRB 051221A. 
However, this puzzle could be resolved if in other directions 
a large amount of baryons (~ 0.01 Mq) ejected in the dou- 
ble neutron star merger have existed there, as found in the 
previous numerical simulations (Rosswog et al. 1999; Ruffert 
& Janka 2001). The ejected material would be accelerated 
by the magnetar wind to a bulk Lorentz factor ~ a few, 
provided that Ei^j ~ 10 52 erg. This wide but only mild- 
relativistic outflow will give rise to a very late (~ 10 6 — 10 7 
s after the burst) multi-wavelength re-brightening. However, 
for typical parameters taken in this Letter (the isotropic en- 
ergy is 1.5 x 10 52 erg and the initial Lorentz factor is 5.0, 
i.e., assuming the material ejected from the merger is about 
1.5 x 10~ 3 Mq), the flux is not bright enough to be de- 
tectable. The emission peaks at t ~ 3 x 10 6 s. The 0.3-10 
keV flux is ~ 1 x 10~ 15 erg s^ 1 cm~ 2 which is marginal for 
the detection of Chandra, and the 8.46 GHz flux is ~ 0.01 
mjy. Both are consistent with the extrapolation of current 
observations (Soderberg et al. 2006; Burrows et al. 2006). 
If an event like GRB 051221A takes place much closer, for 

1 The numerical coefficient 3/4 is adopted because in the ideal 
MHD shock jump condition, the random Lorentz factor of the 
post-shock protons is proportional to 3(r rs h — l)/4 rather than 
Trsh — 1 ( scc eq. (17) of Fan, Wei & Zhang 2004b; see also Fan, 
Wei & Wang 2004a and Zhang & Kobayashi 2005 for numerical 
verification). 
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example, at z ~ 0.1, such very late multi-wavelength re- 
brightening may be detectable for Swift XRT or Chandra 
and other radio telescopes. This predication could be tested 
in the coming months or years. 



4 DISCUSSION AND SUMMARY 

Short hard GRBs may be powered by the merger of double 
neutron stars (e.g. Eichler et al. 1989). For GRB 051221 A, 
ruling out of bright supernova component in the late opti- 
cal afterglow lightcurve suggests that the progenitor prob- 
ably is not a massive star and instead is consistent with 
the double neutron star merger model. After the energetic 
merger, a black hole (Eichler et al. 1989), or a differentially 
rotating neutron star (Kluzniak & Ruderman 1998; Dai et 
al. 2006), or a magnetar may be formed. Among the above 
cases the last one is based on various dynamo mechanisms 
(Rossowog, Ramirez-ruiz & Davis 2003; Gao & Fan 2006 
and references therein) and in particular the MHD simula- 
tion of the two neutron star coalescence (Price & Rosswog 
2006). How long can the supermassive magnetar survive? In 
the general-relativistic numerical simulation, the resulting 
hypermassive magnetar collapses in a very short time ~ 100 
ms (Shibata et al. 2006) . The hypermassive magnetar has a 
mass exceeding the mass limit for uniform rotating but the 
supermassive magnetar does not. Therefore a supermassive 
magnetar may be able to survive in a long time as shown in 
Duez et al. (2006). In view of the uncertainties involved in 
the numerical simulation, some constraints from the obser- 
vation rather than just from the theoretical calculation are 
needed. 

Thanks to the successful running of Swift and Chandra, 
now short GRBs could be localized rapidly and their X- 
ray afterglows could be monitored continuously. The X-ray 
fiat segment in GRB 051221A strongly suggests a significant 
energy injection. Such an energy injection could be well ac- 
counted for if the central engine is a millisecond pulsar with 
a dipole magnetic field ~ 10 14 Gauss. The X-ray flat seg- 
ment in GRB 051221A thus provides us a possible evidence 
for a long time living magnetar formed in the double neutron 
star merger. In this scenario, the material ejected from the 
merger would be swept and accelerated by the strong magne- 
tar wind. This wide but mild-relativistic component would 
give rise to a very late multi-wavelength re-brightening and 
might be detectable for an event like GRB 051221A but 
much closer. 

We would like to point out that the magnetar energy 
injection model is not unique to account for the data. For 
example, assuming the energy carried by the material of the 
initial GRB ejecta satisfies the relation E(> F) oc r~ 4,5 
(Rees & Meszaros 1998; Sari & Meszaros 2000), the X-ray 
afterglow lightcurve of GRB 05 1221 A could also be well re- 
produced (Soderberg et al. 2006; Burrows et al. 2006). How- 
ever, the physical process pumping such kind of energy in- 
jection, in particular in the short GRB scenario, is not clear 
yet. 
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